Spatiotemporal structure of Lyapunov vectors in chaotic coupled-map lattices 
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The spatiotemporal dynamics of Lyapunov vectors (LVs) in spatially extended chaotic systems is studied by 
means of coupled-map lattices. We determine intrinsic length scales and spatiotemporal correlations of LVs 
corresponding to the leading unstable directions by translating the problem to the language of scale-invariant 
growing surfaces. We find that the so-called 'characteristic' LVs exhibit spatial localization, strong clustering 
around given spatiotemporal loci, and remarkable dynamic scaling properties of the corresponding surfaces. In 
contrast, the commonly used backward LVs (obtained through Gram-Schmidt orthogonalization) spread all over 
the system and do not exhibit dynamic scaling due to artifacts in the dynamical correlations by construction. 
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Lyapunov exponents (LEs) measure the exponential sep- 
aration (or convergence) of nearby trajectories and provide 
the most common tool to characterize spatiotemporal chaos 
(STC) |1, 2, 3]. Not only exponential separation rates but 
also spatial correlations are crucial to deal with predictability 
questions in extended systems |4]. Basically, (almost) any ini- 
tial infinitesimal perturbation evolves in time asymptotically 
aligning along the most unstable direction, i.e. the main Lya- 
punov vector (LV). In extended systems, the main LV is found 
to be localized in space at any given time. However, due to 
spatial homogeneity, all directions in tangent space are actu- 
ally equivalent and the localization center keeps moving all 
over the system. This is known as 'dynamic localization' of 
the mainLV Sll]. 

In the last decade, there has been some progress in the 
study of STC with tools borrowed from nonequilibrium sta- 
tistical physics. In particular, the evolution of perturbations in 
spatially extended chaotic systems has been shown to be de- 
scribed by multiplicative Langevin-type equations |l5l|6l|7l|8t]. 
In many cases |6], the dynamics of perturbations can be ex- 
pressed in terms of the prototypical stochastic surface growth 
equation of Kardar-Parisi-Zhang (KPZ), dth = (dxh)^ + 
dxxh + where ^ is a white noise [£]• 

In this Rapid Communication we study the spatiotemporal 
structure of STC encoded by the LVs. We show that a fam- 
ily of vectors — that we shall call characteristic LVs — caiTy 
important information about the real-space structure, localiza- 
tion properties and space-time correlations. These properties 
are disclosed only after a logarithmic transformation, so that 
each LV is mapped into a rough surface. Our results greatly 
strengthen the link between STC and certain nonequilibrium 
surface growth models. We also find that the most widely used 
orthogonal LVs — that appear in the standard Gram-Schmidt 
procedure to compute the LEs ifTol [Till — lack many of these 
features due to the construction procedure and, therefore, have 
much less physical significance. 

We illustrate our results by numerical simulations of cou- 
pled map lattices (CMLs), which are a prototype of spatially 
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extended dynamical systems exhibiting chaos. The main ad- 
vantage of CMLs, as compared with e.g., partial differential 
equations, is their lower computational cost. We consider an 
array of L maps diffusively coupled to nearest-neighbors: 

u^{t + 1) = ef{u,+i{t)) + ef{u,^i{t)) + (1 - 2e)/(u,(i)), 

(1) 

where e is the coupling parameter and / a map with chaotic 
dynamics. Throughout this paper, we consider the logistic 
map f{y) — 4j/(l — y) with periodic boundary conditions. 
We follow the standard convention and sort LEs in decreasing 
order Ai > A2 > • • ■ > Al. The coupling parameter is chosen 
to be e = 0.1, and the system is hyperchaotic with A„ > for 
n/L < 0.795. 

Lyapunov vectors.- Given the initial state of the sys- 
tem u(< = 0) = [ui(0), 1*2(0), ul(0)], any infinitesi- 
mally small change Su{t — 0) in the initial condition evolves 
up to linear order (i.e. in tangent space) according to 

Su,{t + 1) = ef{u.,+i{t))5u,+i{t)+ef{u,^i{t))Su,^i{t) 
+ (1 - 2e)f'{u,it))du,{t) = T,j[u{t)]6uj{t) (2) 

with /' being the derivative of f{y) with respect to y and 
T[M(t)] the L X L Jacobian matrix evaluated at u{t). The 
evolution of an infinitesimal perturbation du{ti) is gov- 
erned by the linear equation (tangent space): ^11(^2) — 
M(i2, ti)Su{ti). The linear operator M is just the product 
of the Jacobian matrices evaluated along the system trajec- 
tory from ti to t2, i.e. M(t2,ti) = T[u(t2 - l)]T[u{t2 - 
2)]...T[u{h)]. 

According to Oseledec's theorem l\jS\ (details can be 
found in (ill) there exists the limit operator $00(^2) = 

limti^_oo[M(t2, ii)M(t2, ti)*] ^('2-'i) , such that the loga- 
rithms of the eigenvalues are the LEs A„, and the eigenvec- 
tors form an orthonormal basis je„(t)|. This set of eigen- 
vectors, so-called backwardLWs 113], indicates the directions 
of growth of perturbations grown since the remote past. The 
backward LVs are precisely the orthonormal vectors obtained 
using the standard Gram-Schmidt orthogonalization method 
to compute the LEs flill . 

As we will show below a much more interesting set 
of vectors {g„(i)} are the so-called characteristic LVs. 
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FIG. 1 : (Color online) (a,b) Snapshot of the profiles of LVs with 
n — 1,2 for system (TJ with L = 512 (note that we are taking the 
absolute values and a logarithmic scale). (c,d) Space-time evolution 
of the maxima of LVs from the 1st to the 10th. The position of the 
maximum for the first LV is indicated with big • symbols. Other 
symbols correspond to LVs for n = 2, . . . , 10. 



These are free-evolving perturbations that grow with expo- 
nent A„ (— A„) when integrating to the far future (past): 



lim 
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,(i2 - h)-^\n\\M{t2,h)g„iti 



A„. Charac- 



teristic LVs are independent of the scalar product and do not 
form an orthogonal basis, in contrast to backward LVs. Also, 
they exhibit remarkable physical properties as they reduce to 
the Floquet eigenvectors for a periodic orbit flS*] ; and for au- 
tonomous continuous-time systems the null LE is associated 
with a characteristic LV tangent to the trajectory. This con- 
trasts with the counterintuitive arrangement of backward LVs. 
We want to emphasize that characteristic LVs evolve freely, 
whereas backward LVs (other than the first one) do not map 
into themselves under time evolution. 

Recently, characteristic LVs are receiving growing atten- 
tion ifisi [isi [TgI [TtIi . Unfortunately, these vectors are difficult 
to compute because information from both the remote past 
(via backward LVs), and the remote future (via forward LVs) 
is needed |21]. We have computed the characteristic LVs from 
the intersection of embedded subspaces from backward and 
forward LVs as described in Refs. ifTsi [l6ll . The nth charac- 
teristic LV is a linear combination of the backward LVs from 1 
ton, and viceversa; and in particular gi (i) = ei(t). Our com- 
putational resources allowed us to achieve the very reasonable 
system size of L — 1024. 

Localization of Lyapunov vectors. - Figure [TJa,b) shows the 
first and second LVs in logarithmic scale; the plots indicate 
strong localization. In Fig. [nc,d) we plot the time evolution 
for the maxima of |e„(t)| and |gn(^)| for n — 1, . . . , 10 back- 
ward and characteristic LVs, respectively. We see that max- 
ima are spread all over the system in the case of backward 
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FIG. 2: (Color online) (a) Spatiotemporal plot in gray scale of \dxr]2 \ 
for the backward LVs. The plateaus of 772 appear as white regions, 
(b) Snapshots of r7„, n = 2, 5, 20 from top to bottom (the curves 
are shifted to improve the visibility), (c) Probability distribution of 
plateau lengths of the same r]n. The dashed straight line is a power 
law ~ . (d) The cut-offs of the plateau distributions as a func- 

tion of n; the (blue) line is a guide to the eye with slope —1.29. 



LVs. However, the localization positions of different charac- 
teristic LVs appear in clusters, so that different vectors tend 
to be highly correlated, indicating that unstable manifolds are 
nearly tangent. This clearly shows that the spread of local- 
ization positions in the case of backward LVs is a byproduct 
of orthogonalization and does not represent properly the spa- 
tial structure of the unstable directions. Extensivity of STC 
has often been interpreted as the independent contribution of 
individual chaotic degrees of freedom located at different po- 
sitions of the system. This has been claimed to be supported 
by the fact that backward LVs do localize at different posi- 
tions flS*]. However, this picture of STC in terms of almost 
independent building blocks, although appealing, becomes 
problematic if we consider the characteristic LVs as the proper 
vectors truly encoding the system dynamics. 

Surface growth picture.- The scaliiig properties of the first 
LV were studied by Pikovsky et al. O, la]. For a wide class 
of extended dynamical systems, they concluded that In |ei it) \ 
can be interpreted as a growing surface exhibiting dynamic 
scaling properties in the universality class of KPZ. In the same 
spirit, here we define a surface for every (backward or charac- 
teristic) LV, f„(t) = [fn{x,t)]%Zi, via the Hopf-Cole trans- 
formation, /i„ (a;, t) = In |/„(a;, Under this logarithmic 
transformation the nth LE maps into the average velocity of 
the corresponding nth surface. 

Some important features of the spatial structure are revealed 
by introducing the auxiliary fields, ?7„ = /i„ — hi, measur- 
ing the relative growth of the nth LV with respect to the first 
one. In Fig. |2lb) we plot snapshots of rjn{x, t) for backward 
LVs with n = 2, 5 and 20 (similar profiles are observed for 
characteristic LVs). The profiles exhibit a distinctive structure 
organized in flat regions. These plateaus evolve so that their 
boundaries have an erratic motion as can be seen in Fig.|2a). 
This spatial structure results from the fact that the dynamics of 
every LV is governed by Eq. (|2|i with the same 'deterministic 
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FIG. 3: (Color online) (a,b) Structure factors of n = 
1, 5, 10, 20, 40 from top to bottom, L = 1024. (c,d) Data collapse 
through scaling relation {3} with fitting parameter 9 = 0.9 ± 0.1. 

noise'. Spatial regions where dxrjnix, i) « correspond to 
regions where the ?ith LV surface closely follows the first LV. 
This means that the dynamics specific of the nth LV surface 
is basically associated with that of the boundaries between 
plateaus. Since the first LV surface dynamics belongs to the 
KPZ universality class, the picture that immediately emerges 
from Fig.|2b) is that the nth LV surface is, loosely speaking, 
"piecewise KPZ". 

The probability distributions of plateau lengths decay in 
good approximation as a power law P(^) oc / for I < A„ 
with 7 « 1.28 ± 0.10 [Fig. I^c)]. We find good evidence 
that the cut-off scales as A„ ~ [L/{n — 1/2)]'^ where x = 
1.29 ± 0.2 [Fig.|2|d)]. There is some numerical uncertainty in 
the plateau and cut-off determination but we can estimate the 
average plateau length (In) « [L/(n - l/2)]0 '^5±"-2. 

Spatial structure.- Let us now study the spatial correla- 
tions of the LV surfaces hn{x, t). We computed the stationary 
structure factor Sn{k) = limt^oo{hn{k, t)hn{-'k, t)), where 
hn{k,t) = '^^exp{ikx)hn{x,t). One expects the first LV 
surface correlations to decay as [Fig. [3|a,b)], in accor- 
dance with KPZ universality f3, Interestingly, the nth 
LV surface for n > 1 also shows scale-invariant correlations 
~ k^^, but only up to a crossover length scale that depends 
on n. At long wavelengths correlations of LV surfaces, as- 
sociated with backward and characteristic LVs, decay as k^^ 
and fc^^'^, respectively (Fig. [3]). This l/fc-divergence indi- 
cates extremely weak long-range spatial correlations for both 
classes of LVs. 

Numerical data for all LV surface spatial correlations in the 
stationary state [Fig. [3ja,b)] can be cast in a single scaling 
function as shown in Fig.[3tc,d): 

Sn{k)L-^k^ = g{k/K) (3) 

where g{u) = const for u 3> 1, and g(u) ^ u"' for 
M <C 1, cr = 1 and a — 0.8 for backward and character- 
istic LVs, respectively. The crossover wavelength is Ik„ ^ 
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FIG. 4: Time evolution of (j)io(k, t) for (a) backward and (b) charac- 
teristic LVs (L = 512). 

[L/{n - 1/2)]-®, with 6* = 0.9 ± 0.1, and gives a typical 
length scale oc for the nth LV surface, which immedi- 
ately suggests to link with the average plateau length (Z„) 
introduced above. 

Space-time correlations.- For each n we consider the field 
hn{x,to + t) — hn{x,tQ), measuring the relative growth of 
the nth LV between times to to to + ^ at ^ given point in 
space and whose structure factor we denote by 4>n{k, t). We 
run trajectories for a long time to to allow a good projec- 
tion into the attractor and then measure (/)„(fc,t) at different 
times t. Figure Ua) shows that backward LV surfaces for 
n > 1 show long-range correlations ^ 1/k already at very 
short times. We claim this is due to non-local information in- 
troduced by the Gram-Schmidt (orthogonalization) procedure 
that involves all LVs with n' < n. Contrastingly, characteris- 
tic LV surfaces [Fig. SJb)] exhibit excellent dynamic scaling 
properties. In this case the nth surface is uncorrected — as in- 
dicated by a flat ~ k'^ structure factor — at scales k ^ t~^/^, 
where z is the dynamic exponent. As time proceeds correla- 
tions develop and extend over a correlation length ^ ^ t^/^, 
as expected for typical dynamic scaling in kinetic surface 
roughening. Eventually, saturation occurs when ^ reaches 
the system size L, and we recover the stationary correlations 

(j)nik, t ^ OO) = Sn{k). 

In agreement with the results for the static correlations, 
Sn{k), discussed above, we can see two scaling regimes in 
Fig. Ub) separated by the typical wavenumber Ik„. Curves 
corresponding to different times in Fig.Ub) can be 'collapsed' 
into: 

</.„(fc,t) = fc-'^'c«(fcti/^') (4) 

where the label i G {a, b} denotes the scaling region, above 
or below the crossover. On the one hand, for t ^ k^^S we 
find Ha = 1.2 [Fig.|5ja)]. Our simulations indicate that the 
dynamic exponent = 1 is consistent with numerics on this 
range of scales for all n. This suggests a ballistic motion of 
the plateau boundaries of This small value of Za implies 
that the surface height correlates very rapidly once the corre- 
lation length reaches the intrinsic length scale Thus, the 
larger n the sooner 4>nik, t) reaches the asymptotic value. On 
the other hand, for t < k"'''' we find fib = '2 [Fig.|3b)]. Our 
simulations indicate — supported also by the computation of 
the surface growth exponent (3^ (not shown) — a slow contin- 
uous increase of the dynamic exponent with n (for n = 1, 
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FIG. 5: (Color online) Dynamic scaling through scaling relation ©. 
(a) Data collapse above the crossover with fia ~ 1-2 and Za — 1. 
Scaling shows up only for wavenumbers smaller than the intrinsic 
scale k < kio ~ 0.13, in agreement with the crossover picture, (b) 
Collapse below the crossover, n — 5 (instead of n — 10) has been 
used to obtain a larger region of scaling; /ij, = 2and2i, = 1.75 ±0.1. 

zkpz = 3/2 for sufficiently large systems). Therefore, at 
short scales characteristic LV surface correlations propagate 
slightly more slowly when we look at larger n. 

Observations and Outlook.- In our opinion characteristic 
LVs are more fundamental for a theoretical understanding of 
STC than the more commonly used (Gram-Schmidt) back- 
ward LVs. We have shown that they have strong localization 
properties, so that they generically localize in physical space, 
signaling regions where most unstable expanding directions 
concentrate at a given time. This suggests that characteristic 
LVs, rather than backward vectors, may play an important role 
in the hierarchical decomposition of STC iflSll . 

We conjecture that the dynamics of the nth characteristic 
LV surface for n > 1 is described by a KPZ saddle solution: 
given a certain realization of the (deterministic) noise there 
should exist (measure zero) initial perturbations that evolve 
toward states whose space-time surface coiTelations scale with 
critical exponents that are not given by the standard KPZ uni- 
versality class. This would coiTespond to the scaling ^ fc^^^ 
for long wavelengths, k < k„, in Fig. [3b). While solutions 
whose basin of attraction is a set of measure zero are irrele- 
vant for stochastic equations, in deterministic chaos these may 



be fundamental solutions for the complete understanding of 
sensitive dependence in high-dimensional chaotic systems. In 
fact, it has recently been shown that the stochastic KPZ equa- 
tion does exhibit this type of saddle solutions for noise realiza- 
tions with suitable (and admittedly bizarre) coiTelations ifioll. 
We have shown that, in statistical terms, the nth characteris- 
tic LV surface hn{x, t) differs from a genuine KPZ surface in 
an amount given by a field rin{x, t). It is not difficult to con- 
vince oneself that the step field corresponds to the nontrivial 
solutions {i.e. rjn const.) of the stochastic equation: 

= {dxVn)'^ + da^xT^n + 2 (Sj: ft-i ) ((9a;?7n ) ■ (5) 

In a quite natural way the term dxrin on the right hand side 
leads to solutions of the type shown in Fig. |2|b); plateaus of 
random spatial extent separated by boundaries moving eiTat- 
ically. We believe that Eq. ^ may indeed provide a mathe- 
matical tool to address questions like real-space structure of 
the dynamics of extended systems, space-time coiTelations in 
the propagation of disturbances, degree of localization of un- 
stable directions, and extensivity and decomposition of STC. 

Finally, we have found identical behavior in other CML 
models. In fact, our conclusions should be generic for LVs 
corresponding to LEs well above zero in systems whose 
first LV surface belongs to the KPZ universality class f20t]. 
(This certainly should include all those systems reported 
in Ref. |@].) Although more research is needed to test the 
generality of our findings, in particular for continuous-time 
systems, our preliminary results with a minimal stochastic 
model (dtw = (w + dxxw) confirm the universality of our 
results l,20il . 
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